Under consideration for publication in J. Fluid Mech. 

The Universal Aspect Ratio of Vortices in 
Rotating Stratified Flows: Theory and 

^ ■ Simulation 

o. 

^' PEDRAM HASSANZADERi, PHILIP S. MARCUSH 

C : AND PATRICE LE GAL2 






in 

rn 



X 



^Department of Mechanical Engineering, University of California, Berkeley, CA 94720, USA 

^Institut de Recherche sur les Phenomenes Hors Equilibre, UMR 7342, CNRS - Aix-Marseille 
Universite, 49 rue F. Joliot Curie, 13384 Marseille, Cedex 13, France 



^ ' (Received ?; revised ?; accepted ?. - To be entered by editorial office) 

>>; 

' O , We derive a relationship for the vortex aspect ratio a (vertical half-thickness over hori- 

^ ■ zontal length scale) for steady and slowly evolving vortices in rotating stratified fluids, as 

;^ , a function of the Brunt- Vaisala frequencies within the vortex Nc and in the background 

fluid outside the vortex TV, the Coriolis parameter /, and the Rossby number Ro of the 
Q ', vortex: a^ = Ro{l + Ro) p/{N^ — N^). This relation is valid for cyclones and anticyclones 

^ ' in either the cyclostrophic or geostrophic regimes; it works with vortices in Boussinesq 

J>-,, fluids or ideal gases, and the background density gradient need not be uniform. Our re- 

rS^ • lation for a has many consequences for equiUbrium vortices in rotating stratified flows. 

, Mh | For example, cyclones must have N^ > N"^; weak anticyclones (with \Ro\ < 1) must have 

A*"^ < 7V^; and strong anticyclones must have N^ > N'^. We verify our relation for a with 
04 \ numerical simulations of the three-dimensional Boussinesq equations for a wide variety 

^ ' of vortices, including: vortices that are initially in (dissipationless) equilibrium and then 

evolve due to an imposed weak viscous dissipation or density radiation; anticyclones cre- 
ated by the geostrophic adjustment of a patch of locally mixed density; cyclones created 
by fluid suction from a small localised region; vortices created from the remnants of the 
violent breakups of columnar vortices; and weakly non-axisymmetric vortices. The values 
f— ■^ \ of the aspect ratios of our numerically-computed vortices validate our relationship for 

^sj ' a, and generally they differ significantly from the values obtained from the much-cited 

conjecture that a = f / N in quasi-geostrophic vortices. 
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1. Introduction 

Compact three-dimensional baroclinic vortices are abundant in geo- and astrophysical 
flows. Examples in pl anetary atmospheres include the rows of cyclor ies and anticyclones 
near Saturn's Ribbon (jSavanagi et a/.ll2010l ) and nea r 41°S on Jupi ter (JHumphrevs fc MarcusI 



20071 ). and Jupiter's anti cyclonic Gr eat Red Spot ('Marcus"l993|). In the Atlantic Ocean 



meddles persist for years (jArmi et al. 1988; Mc Williams 1985) , and numerical simulation s 



of the disks around protostars produce compact anticyclones ([Barranco fc Marcusll2005l ). 
The physics that create, control, and decay these vortices is highly diverse, and the as- 
pect ratios a = H/L of these vortices range from flat "pancakes" to nearly round (where 

f Email address for correspondence: pmarcus@me.berkeley.edu 
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H is the vertical half- height and L is the horizontal length scale of the vortex). However, 
we shall show that the aspect ratios of the vortices all obey a universal relationship. 
Our relation for a differs from previously published ones, including the often-used 

a — f /N, where / is the Coriolis parameter, N = y— ^ gf is the Brunt- Vaisala frequency, 

g is the acceleration of gravity, z is the vertical coordinate, p is the density for Boussinesq 
flows and potential density for compressible flows; and a bar over a quantity indicates 
that it is the value of the unperturbed (i.e., with no vortices) background flow. We shall 
show that a = f / N is not only incorrect by factors of 10 or more in some cases, but also 
that it is misleading; it suggests that a depends only on the background flow and not 
on the properties of the vortex, so that all vortices embedded in the same flow (e.g., in 
the Atlantic or in the Jovian atmosphere) have the same a. We shall show that this is 
not true. Knowledge of the correct relation for a is important. For example, there has 
been debate over whether the color ch ange, from white to r ed, of Jupiter's anticyclone 
Oval BA, was due to a change in its H ( de Pater et a^.lboiof ). Measurements of the half- 



heights H of planetary vortices are difficult, but H can be accurately inferred from the 
correct relation for a. We validate our relation fo r a with 3D numeri cal simulations of 
the Boussinesq equations. A companion paper bv lAubert et al\ ( 20121 ) validates it with 



laboratory experiments and with observations of Atlantic ocean meddles and Jovian 
vortices. 



2. Aspect Ratio: Derivation 

We assume that the rotation axis and gravity are parallel and anti-parallel to the 
vertical z axis, respectively. We also assume that the vortices are in approximate cyclo- 
geostrophic balance horizontally and hydrostatic balance vertically (referred to hereafter 
as CG-H balance). Necessary approximations for CG-H balance are that the vertical Vz 
and radial Vr velocities are negligible compared to the azimuthal one vg (where the origin 
of the cylindrical coordinate system is at the vortex centre), that dissipation is negligible, 
and that the flow is approximately steady in time. With these approximations, the radial 
r and azimuthal 9 components of Euler's equation in a rotating frame are 

dp/dr — pvg{f + vg/r) and dp/dz = —pg, (2.1) 

where p is the pressure. We have assumed that the vortex is axisymmetric, but will show 
later numerically that this approximation can be relaxed. Following the convention, we 
ignored the centrifugal term p/^r/4 f in equation ()2.ip by assuming that the centrifu- 
gal buoyancy is much smaller than the gravitational buoyancy, i.e. that the rotational 
Froude number pd/{4:g) <C 1, where d is the charac teristic distance of the vortex from 
the rotation axis (see e.g. iBarcilon fc Pedloskvlll967 ). The 0-component of Euler's equa- 



tion, continuity equation, and the equation governing the dissipationless transport of 
(potential) density are all satisfied by a steady, axisymmetric flow with Vr — Vz — 0. As 
a consequence, equations (|2.1|) are the only equations that need to be satisfied for both 
Boussinesq and compressible flows. Thus, our relation for a will also be valid for both of 
these flows. Far from the vortex, where v = 0, p = p, and p — p, (|2.ip reduces to 

dp/dr = and dp/dz = —pg (2-2) 

showing that p and p are only functions of z. Subtracting equations (|2.2p from (12. 1[) : 

dp/dr = pvg (/ + vg/r) (2.3) 

dp/dz = -pg (2.4) 
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where p = p — p and p = p — p arc respectively the pressure and density anomahes. The 
centre of a vortex (r = z = 0) is defined as the location on the z-axis where p has its 
extrenium, so equation (12.41) shows that at the vortex centre (denoted by a c subscript) 
Pc = OT pc = p(0) = Po- At the vortex boundary and outside the vortex, where ve and 
p are negligible, equations (j2.3p and (|2.4p show that p ~ 0. 

We define the pressure anomaly's characteristic horizontal length scale (i.e. radius) 
as L = a/|4pc/(V^p)c|, where the subscript _L means horizontal component. Integrating 
(|2.3p from the vortex centre to its side boundary at (r, z) = {L, 0) approximately yields 

-Pc/L^PoVeif + Ve/R,) (2.5) 

where in the course of integration, p has been replaced with po, which is exact for Boussi- 
nesq flows, and an approximation for fully compressible flows. Here Ve is the characteristic 
peak azimuthal velocity, and Rv is the approximate radius where the velocity has that 
peak. The analytical and numericall y simulated vo r tices discussed below, meddles, and 
the laboratory vortices examined bv lAubert et ali ( 20121 ) all have Ry — L, but hollow 



vortices wit h quiescent interiors have Ry ^ L. For example, the Great Red Spot has 
Ry ~ 3i (Shettv fc MarcuJbOld ). Similarly, integrating (|2.4p from the vortex centre to 
its top boundary near (r, z) = (0,7J) approximately gives 

-Pc/H^-gp{r^O,z^H), (2.6) 

where H is the pressure anomaly's characteristic vertical length scale (i.e. half- height), 
H=^\2p,/{d4/dz^),\. 

Equations (|2.5p and (|2.6p can be combined to eliminate Pc'- 



PoVeif + Ve/Rv) ^ gp(r = Q,z^H) 

H L ' ^ ' ' 

Notice that this equation is basically the ther mal wind equ ation, with the cyclostrophic 
term included (i.e. the gradient- wind equation (|Vallisll2006[ )'). integrated over the vortex. 



Using the first term of a Taylor series, we approximate p{r = 0, z = H) on the right-hand 
side of (l2Jl) with 



p{r = 0,z = H) = p, + H{dp/dz), = H[{dp/dz), - [dp/dz),] = poH{N^ - N^)/g (2.8) 

where p = p — p and pc = have been used. Note that in general, N{z) is a function of z; 
however, the only way in which N{z) is used in this derivation (or anywhere else in this 
paper) is at z = for evaluating {dp/dz)c- Therefore, rather than using the cumbersome 
notation Nc, we simply use N. 

Using (|2.8p in equation (|2.7p gives our relation for a: 

2_fHY Ro[l + Ro{L/Ry)] 2 



\L ) N^-N^ -^ ^ ' 

where the Rossby number defined as Ro = Ve/{fL) can be well approximated as Ro = 
Wc/(2/), ijJc being the vertical component of vorticity at the vortex centre. Defining the 
Burger number as Bu = (N H / {f L))'^ , equation (|2.9p may be as well rewritten as 

[{NJN)^ - l]Bu = Ro{l + {L/Ry)Ro) (2.10) 

Equation (J2.9I) shows that a depends on two properties of the vortex: Ro and the 
difference between the Brunt- Vaisala frequencies inside the vortex (i.e. N"^) and outside 
the vortex (i.e. N"^). Note that to derive relation (j2.9p . no assumption has been made 
on the compressibility of the flow, Rossby number smallness, dependence of N on z, or 
the magnitude of Nc/N. Therefore, equation (|2.9p is applicable to Boussinesq, anelastic 
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( VallisI 120061 ) ■ and fully compressible flows, cyclones (i.e., Ro > 0) and anticyclones (i.e., 



Ro < 0), and geostrophic and cyclostrophic flows. In the cyclostrophic limit (i.e., \Ro\ ^ 
1) with Nc = and Ry = L Ro(l+Ro) -^ i?o^ henc e equation (12.91) becomes Ve = HN, 
agreeing with the findings of lBillant fc Chomaa (|200lh and others. Equation (|2.9p is easily 



modified for use with discrete layers of fluid rat her than a conti nuous stratific ation, and 
in that case agrees with the theoretical work of lNoJ (|l98l[) and ICartonI (J2001I ). 

Equation (|2.9I) has several consequences for equilibrium vortices. For example, because 
the right-hand side of (|2.9p must be positive, cyclones must have N^ ^ N'^. Another 
consequence is that anticyclones with ~Ro < Ry/L, must have N^ ^ N'^, and anti- 
cyclones with —Ro > Rv/L, have N^ ^ iV^. In addition, equation (|2.9p is useful for 
astrophysical and geophysical observations of vortices in which some of the vortex prop- 
erties ar e difficult to measu re. For example, Nc is difficult to measure in some ocean 
vortices ( Aubert et aLl 120121 ). and H is difficult t o determine in some satellite observa- 
tions of atmospheric vortices (|de Pater et aLll2010i) . but their values can be inferred from 
equation (12.91) . 

Note that Nc is a measure of the mixing within the vortex; if the density is not 
mixed with respect to the background flow, then Nc ^ N (and the vortex is a tall, 
barotropic Taylor column) ; if the density is well- mixed within the vortex so t he (potential) 



densit y is uniform inside the vortex, then Nc — >■ (as in the experiments of lAubert et i 
( 20121 )): if Nc > N^ (as required by cyclones), then the vortex is more stratified than 
the background flow. 



3. Previously Proposed Scaling Laws 

Other relations for a that differ from our equation (j2.9p have been published previ- 
ously, and the most frequently cited one is a = H/L = f / N. This relationship is inferred 
fr om Charney's e quation for the quasi-geostrophic (QG) potential vorticity (equation (8) 
in ICharnevlll97ll ) that was derived for flows with \Ro\ <^ 1 and Nc/N ~ 1. Separately 
re-scaling the vertical and horizontal coordinates of the potential vorticity equation, 
and then assuming that the the vortices are isotropic in the re-scaled (but not physi- 
cal) coordinates, one obtains the alternative scaling a — f/N. Numerical simulations 
of the QG equatio n for some initial conditions have produced turbulent vortices wit h 
H/L « f/N (c.f., iMcWilliams et oi.lll999t [Pritschel et qIII1999I: iReinaud et olll2003l ). 
even though significa nt anisotropy in the re-scaled coordinates was observed in similar 
simulations (McWilli ams et al\\l994) . The constraints under which the QG equation is 
de rived are very restric tive; for example, none of meddles or laboratory vortices studied 
bv I Aubert et all ( 20121) meet these requirements because Nc/N is far from unity. There- 
fore, it is not surprising that none of these vortices, including the labora tory vortices, 
agree with a « f/N, but instead have a in accord with relation p.9p (JAubert et al 
2012[) . 



The constraints under which our equation (|2.9p for a is derived are far less restric- 
tive than those used in deriving Charney's QG equation (and we never need to assume 
isotropy). In particular, one of several constraints needed for deriving Charney's QG 
equation is the scaling required for the potential temperature (his equation (3)), which 
written in terms of the potential density is 



p/p=-{f/9)idi'/d^) 



(3.1) 



where ip is the stream function of horizontal velocity. This constraint alone (which is 
effectively the thermal wind equation) implies our relationship ()2.9|) for a. To see this, 
in equation p.ip replace i{j with VgL, d/dz with 1/iJ, Vg with RofL, and p/p with 
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H{{dp/dz)c ~ [dp/dz)}/p = {H/g){N'^ - N^). With these replacements, equation dHU 
immediately gives 

{H/Lr = Rof/{N!-N'), (3.2) 

w hich i s the small Ro limit of equation (|2.9p . 



Gilll ()1981l) also proposed a relationship for a that differs from ours. He based his 
relation for a on a model 2D zonal flow (that is, not an axisymmetric vortex, but rather 
a 2D vortex) and found that a was proportional to Ro f / N . To determine a, Gill derived 
separate solutions for the flow inside and outside his 2D model vortex, which he assumed 
was dissipationless and in geostrophic and hydrostatic balance. Despite the fact that 
Gill's published relation for a, obtained from the outside solution, differs from ours, we 
can show that his solution for the flow inside his 2D vortex satisfies our scaling relation 
for a. Gill's solution for the zonal velocity (which is in the y direction) is u = —{f/a)x 
(his equation (5.14) in dimensional form). His density anomaly is p = po{N'^/g)z (i.e. 
within the 2D vortex, p — po). The equation for v gives oj^ = —{f/a), and therefore 
Ro = uJc/(2f) = — l/(2a). Substituting v and p into the equations for geostrophic and 
hydrostatic balance, gives dp/dx — —po{f^/a)x and dp/dz = —poN^, respectively. 
Using the definitions of H and L from section [2] along with Ro = — l/(2a), we obtain 
a^ = —Rof^/N^, which is our relation (|2.9p in the limit of small Ro, L — Ry, and Nc = 
(which are the constraints under which Gill's solution is obtained). Gill's scaling for a 
is derived from the fiow outside the vortex, which he derived by requiring that both the 
tangential velocity and density are continuous at the interface between the inside and 
outside solutions. In general, this over- constrains the dissipationless fiow (which only 
requires pr essure and no r mal v elocity to be continuous) - see for example the vortex 



solution in lAubert et al\ (|2012| ) in which the pressure and normal component of the 
velocity are continuous at the interface, but not the density or tangential velocity. The 
extra constraints force the solution outside Gill's vortex to have ad ditional, (unphysical) 



length scales, resulting in Gill's relation for a differing from ours. lAubert et al\ ( 20121) 



show that Gill's relationship for a does not fit their laboratory experiments, meddles, or 
Jovian vortices. We examine the accuracy of both Charney's and Gill's relationships in 
section [T] 



4. Gaussian Solution to the Dissipationless Boussinesq Equations 

It is easy to find close d-form solutio n s to t he steady, axisymmetric, dissipationless 
Boussinesq equations fe.g. lAubert et al\ (|2012f) ). One solution that we shaU use to gen- 



erate initial conditions for our initial-value codes is the Gaussian vortex with p = 
Pc exp[— (z/i?)^ — {r/Vf'] and Vr = v^ = Q (where Pc, H, and L are arbitrary con- 
stants). Then, p is found from dp/dz using equation (|2.4|) . and vg is found from dp/dr 
using equation (|2.3|) with p replaced by po- This Gaussian vortex exactly obeys our rela- 
tionship (|2.9p for a when the Rossby number is defined as before as Ro = a;c/(2/), when 
Rv is set equal to L, and when the vertical and horizontal scales are defined as in sec- 
tion[51 Note that N{r, z) within the vortex is not uniform, that 7V^ — N^ ~ 2pc/{poH'^), 
and that the vortex is shielded. By shielded, we mean that there is a ring of cyclonic 
(anticyclonic) vorticity around the anticyclonic (cyclonic) core in each horizontal plane, 
and therefore at each z, circulation due to the vertical component of the vorticity is zero 
(i.e. the vortices are isolated). The Gaussian vortex could be a cyclone or an anticyclone 
depending on the choice of constants. This vortex is well-st udied and has been widely 
used to m odel isolated vortices , especially in t he oceans (e.g. iGent fc McWilliamslll986l : 



Morel fc McWilhams.,1997. : .Stuart et aLll201lh . 
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5. Numerical Simulation of the Boussinesq Equations 

We have used 3D numerical simulations to verify our relation (|2.9|) for a in a Boussinesq 
flow with constant N and /. We include diss ipation and solve the equations in a rotating- 
frame in Cartesian coordinates ( Vallisll2006l ): 



V.v = 0, ^^-^+vx/z-A,, + ,v^., f =^o^--A (5.f) 

Ut Po Po P>t g Trad 

where D/Dt = d/dt + v ■ V, and v = {u,v,w) (Notice that throughout this paper, we 
use Vz and w for the vertical component of velocity in the cylindrical and Cartesian 
coordinates, respectively.) We include kinematic viscosity v, but neglect the diffusion 
of density because diffusion is slow (e.g., for salt-water the Schmidt number is ^ 700). 
Instead, inspired by astrophysical vortices (e.g., Jovian vortices or vortices of protoplan- 
etary disks) for which thermal radiation is the main dissipating mechanism, we have 
added the damping term —p/rrad to the density equation to model radiative dissipation 
where Trad is radiative dissipation time scale. 

A pseudo-spectral method with 256^ modes is used to solve equations ()5.ip in a triply 
periodic domain (which was chosen to be 10 to 20 times la rger than the vortex i n eac h 
direction). Details of the numerical method is the same as iBarranco fc Marcus! ( 20061 ) . 



The results of our triply periodic code are qualitatively, and in most cases quantitatively, 
the same as solutions we obtained with a code with no-slip vertical boundary conditions. 
That is because our vortices are far from the vertical boundaries, and therefore the Ekman 
circulation is absent. 



6. Numerical Results for Vortex Aspect Ratios 

As shown in table [H we have examined the aspect ratios of vortices in four types of 
initial-value numerical experiments. The goal of these simulations is to determine how 
well the aspect ratios a of vortices obey our relation (|2.9p as they evolve in time. 

6.1. Case A: Run-Down Experiments 

In this case, our initial condition is the velocity and density anomaly of the Gaussian vor- 
tex from section[4]that is an exact equilibrium of the dissipationless Boussinesq equations 
with constant / and N . These are "run-down" experiments because they are carried out 
either with radiative dissipation (i.e., finite Trad) or viscosity, but not both. Due to the 
weak dissipation, the vortices slowly evolve (decay) and do not remain Gaussian. Also, 
as a result of the dissipation (and decay), a weak secondary flow is induced (i.e. non-zero 
Vr and Vz)- 

6.2. Case B: Vortices Generated by Geostrophic Adjustment 
This case is motivated by vortices produced from the geostrophic adju stment of a locally 



mixed patch of de nsity, e.g. generated from diapycnal mixing (see e.g. iMcWillia ms 198 



Stuart et aLll201lh . Our flow is initialised with v = and p ^ 0. For Cases Bl and B3 
the initial p is that of the Gaussian vortex discussed in section SI But here, the initial 
flow is far from equilibrium because v = 0. In Case B2, the initial p is Gaussian in r, 
but has a top-hat function in z (for this case, the initial H is defined as the half-height 
of the top- hat function). It is observed in the numerical simulations that geostrophic 
adjustment quickly produces shielded vortices. 
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Case N/N^ 



H/L 



Ro £'fe(10-5)|Case N/N^ 



H/L 



Trad Ro -Bfc(10 



Al 


1/0 


16/8 


-0.2 


25 


A18 


3.33/3.17 


1.44/0.72 


oo 


-0.22 


4 


A2 


1/0 


16/8 


-0.2 


12.5 


A19 


10/10.5 


0.16/0.08 


CX) 


+0.87 


250 


A3 


1/0 


16/8 


-0.2 


6.25 


A20 


3.33/3.49 


1.44/0.72 


1656 


+0.14 





A4 


1/0.5 


16/8 


-0.14 


25 


A21 


1/1.0475 


16/8 


1656 


+0.015 





A5 


1/0.5 


4/2 


-0.14 


25 


A22 


3.33/3.17 


1.44/0.72 


1237 


-0.22 





A6 


20/19.99 


0.05/0.02 


-0.11 


500 


A23 


3.33/3.49 


1.44/0.72 


120 


+0.14 





A7 


1/0 


20/12 


-0.12 


11.1 


A24 


3.33/3.17 


1.44/0.72 


120 


-0.22 





A8 


1/0 


16/8 


-0.2 


25 


Bl 


1/0 


20/12 


oo 


-0.12 


4 


A9 


1/0 


12/12 


-0.04 


11.1 


B2 


1/0 


20/12 


oo 


-0.12 


4 


AlO 


20/19.9 


0.03/0.03 


-0.2 


222.2 


B3 


1/0 


20/12 


oo 


-0.12 





All 


20/19.9 


0.03/0.03 


-0.2 


11.1 


CI 


1/0.761 


6.09/5.67 


oo 


+0.33 


49.7 


A12 


3.53/2.5 


0.96/0.96 


-0.5 


39.3 


C2 


1/0.709 


5.91/5.67 


oo 


+0.35 


45.9 


A13 


3.54/3.45 


0.96/0.96 


-0.023 


39.3 


Dl 


2.5/2.5 


40/2 


oo 


-0.75 


25.4 


A14 


1.56/1.64 


6.55/3.28 


+0.0386 


3.91 


D2 


2.5/2.5 


40/2 


oo 


-0.75 


12.7 


A15 


0.5/0.75 


16/8 


+0.0477 


4 


D3 


1.67/1.67 


90/4.5 


oo 


-0.5 


16.9 


A16 


0.5/0.55 


16/8 


+0.0083 


4 


D4 


5/5 


10/0.5 


oo 


-1 


50.7 


A17 


3.33/3.49 


1.44/0.72 


+0.15 


4 


D5 


20/20 


0.6/.03 


oo 


+0.8 


24.4 



Table 1. Parameters of the background flows and of the vortices at the "initial" time. For 
Cases A, B, and D the "initial" time is f = 0, and for Case C the "initial" time is f = toff- All 
values are in CGS units. For all cases, / = 5 rad/s, g = 980 m/s , and Po = 1 g/cm . Ekman 
number is defined as Ek = v/{fL^). See text for the difference between Cases Bl and B2. For 
Case CI, toff = 60 s and Q — —64 cm'^/s, and for Case C2, toff = 30 s and Q — —128 era? /s. 



6.3. Case C: Cyclones Produced by Suction 



Injection of fluid into a rotating flow generates anticyclones (jAubert et a/.ll2012l ). while 
suction produces cyclones. We simulate suction by modifying the continuity equation in 
(|5.1|) as V • V = Q(x, t) where Q is a specified suction rate function. The flow is initialised 
with V = p = 0. Suction starts at i = over a spherical region with radius of 6 cm and 
is turned off at time toff- A shielded cyclone is produced and strengthened during the 
suction process. As mentioned at the end of section [2l for Ro > 0, relation (|2.9p requires 
Nc > N, which we have shown in the numerical simulations that the initial suction 
creates. Cases CI and C2 have different suction rates and toff, but the same total sucked 
volume of fluid, and it is observed that the produced cyclones are similar. 

6.4. Case D: Vortices Produced from the Breakup of Tall Barotropic Vortices 
The violent breakup of tall barotropic (^-inde pendent) vortices in rotating , stratified flows 



can produce stable compact vortices (see e.g. lSmvth fc McWilliamslll998[ ). In Case D, our 
flows are initialised with an unstable 2D columnar vortex with vg = Ro f rexp{—{r/L)^) 
and p — (for this case, the initial H is the vertical height of the computational domain). 
Note that the initial columnar vortex is shielded. Noise is added to the initial velocity 
field to hasten instabilities. The vortex breaks up and then the remnants equilibrate to 
one or more compact shielded vortices (in each case, only the vortex with the largest 
\Ro\ is considered in section [7]). 



7. Aspect Ratio: Numerical Simulations 

In all cases, vortices reach quasi-equilibrium and then slowly decay due to viscous or 
radiative dissipation except for Case B3 which is dissipationless and evolves only due to 
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Figure 1. Comparison of qnum with qthr (see text for definitions). The circles show the 
value of |_Ro(l + -Ro)|(aNUM/aTHR)^ and the straight lines show the value of this expression if 
aNUM = aTHR- AH 4122 data points (circles) collapse on the straight lines (and densely cover 
them), validating our equation (|2.9p . Data points are recorded one inertial period (= 47r//) after 
the initial time (as defined in table[T]| in Cases A-C, and 50 inertial periods after t = in Case D. 
Note that all of our simulated vortices have L = Rv The horizontal axis in the inset is the same 
as in the main figure; the inset's vertical axis is the relative difference |1 — (oNuivi/aTHR)^! (which 
is < 0.07). {n.b., the left-most plotted point has Ro{l + Ro) ~ —0.25 due to the mathematical 
tautology that Ro{l + Ro) ^ —0.25 for all values of Ro.) 



geostrophic adjustment. As a result, Ro decreases, and the mixing of density in the vortex 
interior changes (i.e., N^ changes). Therefore, it is not surprising that the aspect ratio a 
also changes in time. Quasi-equilibrium is reached in Case A almost immediately. In Case 
B, vortices quickly form and come to quasi-equilibrium after geostrophic adjustment. 
Quasi-equilibrium is achieved following the geostrophic and hydrostatic adjustments after 
toff in Case C, and (much longer) after the initial instabilities in Case D. 

For each case, we use the results of the numerical simulations to calculate Ro{t) = 
ujc{t)/{2f) and N^it) = y/N^ - g{dp{t)/dz)J p^. We compute L{t) and H{t) from the 
numerical solutions using their definitions given in section [5] Calculating L based on V^ 
rather than just r-derivatives is useful for non-axisymmetric vortices. For example, due 
to a small non-axisymmetric perturbation ad ded to the initial condition of C ase A8, the 
vortex went unstable and produced a tripole ( Van Heiist fc Kloosterziellll989l ). Cases CI 
and C2 also produced non-axisymmetric vortices. We define the numerical aspect ratio 
as aNUM(i) = H{t)/L{t). We define the theoretical aspect ratio aTHR from equation ()2.9p 
using Ro{t) and Nc{t) extracted from the numerical results and the (constant) values of 
/ and N. 

Figure [T] shows how well axHR agrees with aNUM • The inset in figure [1] shows that 
the relative difference between the two values, calculated as |1 — (aNUM/axHR)^!, is 
smaller than 0.07. For each case, the maximum difference occurs at early times or during 
instabilities. 

Figure [2] compares the values of aNUM with axHR as a function of time for six cases. 
The figure starts at time t = 0, so it includes vortices which are not in CG-H equilibrium 
to highlight the situations for which relationship (j2.9p for a is not good due to viola- 
tion of its assumptions. Cases Al, Bl, and A20 in figure [2al exhibit excellent agreement 
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with our theoretical prediction for a, while Case A8 shows a small deviation starting 
around t = 80(47r//). This deviation is a result of the vortex going unstable at this time 
(accompanying by relatively large w,- and Vz) and forming a tripolar vortex . After the 
tripole comes to CG-H equilibrium, its a once again agrees with theory. As the vortices 
dissipate, and Ro and Nc change, a can either decrease in time (c.f.. Case Al) or increase 
(c.f.. Case A20). 

Figure [2bl shows Cases Dl and D3 from time t = 0. The remnant vortices that formed 
from the violent break-up of the columnar vortices are initially far from the CG-H balance. 
As a result, the value of q;thr at these early times does not fit well with the values of 
Q^NUM- However, after the CG-H balance is established in the remnants, our theoretical 
relationship (|2.9p becomes valid and axHR agrees well with qnum- 

Figure [2a] shows that the alternative scaling relation based on Charney's QG equation, 
a — f/N, is not a good fit to our numerical data. Cases Al, A8, and Bl all have 
= 5 which is obviously far from the measured aspect ratio of these vortices. Case A20 
has f/N = 1.5 which again does not agree with onum- In fact, in all four cases, the 
difference between a-^vMit) and f/N increases by time, while aTHR(0 always remains 
close to a-^uuit)- For other cases in table [TJ it has been observed that for vortices which 
are in CG-H equilibrium, aNUM/(//-^) can be as large as 9.56 and as small as 0.11. The 
data displayed in figure I2bl were carefully "cherry-picked" from all of our runs because 
they are unusual in that a — t- f/N after a long time. The fiuid within the remnants 
strongly mixed with the background fiuid, so at late times Nc -^ N and Ro significantly 
decreases and therefore the conditions needed for the validity of Charney's QG equation 
are approached. Whether these results are a fiuke and whether a —^ f/N for all vortices 
that are created by one particular method is not yet clear. The physics governing these 
vortices is currently be investigated and will be discussed in a future paper. 

Gill's model ( jGiH 1981) , discussed in section |3l is not a good fit to any of our nu- 
merically computed vortices. For example, the value of at^uuit) / {Ro{t) f / N) is between 
2 and 8 for Case Al; 2 and 9 for Case A8; 80 and 160 for Case A20; and 4 and 7 for 
Case Bl. The much larger error observed for Case A20 is due to the fact that unlike the 
other three cases, Nc is far from in this case, and Gill's derivation does not incorporate 
Nc^O. 



8. Conclusion 

We have derived a new relationship (|2.9p for the aspect ratio a of baroclinic vortices in 
cyclo-geostrophic and hydrostatic (CG-H) equilibrium and used numerical initial-value 
simulations of the Boussinesq equations to validate this relation for a wide variety of 
unforced quasi-steady vortices generated and dissipated with different mechanisms. Our 
new relationship shows that a depends on the background flow's Coriolis parameter / 
and Brunt- Vaisala frequency N, as well as properties of the vortex, including Ro and Nc- 



same aspect ratios. In a companion paper. lAubert et al\ ()2012l) verify the new relationship 



Kgn 

m 



with laboratory experiments and show it to be consistent with observations of Atlantic 
meddles and Jovian vortices. 

Equation (|2.9p for a has several consequences. For example, it shows that for cyclones 
{Ro > 0), Nc must be greater than N, that is, the fluid within a cyclone must be super- 
stratified with respect to the background stratiflcation. Mixing usually de-stratifles the 
flow over a local region, and therefore cannot produce cyclones. This may explain why 
there are more anticyclones than cyclones observed in nature. We numerically simulated 
local suction to create cyclones, and we found that suction creates a large envelope of 



10 



p. Hassanzadeh, P. S. Marcus, and P. Le Gal 




1000 1200 1400 



Figure 2. Time evolution of othr (continuous curves), and q:num for Cases (a) Al o, A8 D, 
Bl ■, and A20 o, (b) Dl o, and D3 ■. Cases Al, A8 and Bl have f/N = 5 and Case A2_0 has 
f/N = 1.5 which differ significantly from a. Case Dl has f/N — 2 and Case D3 has f/N — 3 
which agree with a only at the late times. 



super-stratified flow around the location of the suction and when the suction is stopped, 
the CG-H adjustment makes cyclones. Details of these simulations and results of an 
ongoing laboratory experiment will be presented in subsequent publications. 

It is widely quoted that vortices obey the quasi-geostrophic scaling law a = f/N (i.e. 
Burger number Bu = 1). This is inconsistent with our relationship which written in 
terms of Bu is Bu — Ro{l + Ro)/[{Nc/N)^ — 1]. We found that, with the exception of 
one family of vortice s, the quasi-geostrop hic scaling law was not obeyed by the vortices 
studied here (and bv lAubert et all (120121) ). a r id cou ld be incorrect by more than a factor 
of 10. Another relationship proposed by ICilll (|l981[ ) was also found to produce very poor 
predictions of aspect ratio. 

We found that a can either increase or decrease as the vortex decays, and our relation- 
ship ()2.9p shows that the dependence of a on Nc is specially sensitive when Nc is at the 
order of N, as it is for meddles and Jovian vortices (JAubert et al\\2012i \. Our simulations 
showed that Nc was determined by the secondary circulations within a vortex and that 
those circulations are controlled by the dissipation. In a future paper we shall report on 
the details of how dissipation determines the secondary flows and the temporal evolution 
of Nc, both of which are importa nt in planetary atmosphere s, oceanic vortices, accretion 
disk flows, and planet formation (JBarranco fc Marcusll2005[ ). 
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